
jgLZ 

(PAGES) 


(NASA CR OR TMX OR AD NUMBER) 


* ' I 


(THRU) 

/ 

^(CODE) 

(CATEGORY) 



X-640-65-363 



TW 

l. 

£ &»£ 


ff 

iC 


5533 0 


DAMPING OF 

INTERSTELLAR PLASMA WAVES 
BY THE COSMIC RAY GAS 


BY 

DEREK A. TIDMAN 


GPO PRICE $. 
CFSTI PRICE(S) S 

Hard copy (HC) 
Microfiche (MF) 

ff 653 July 65 


j 



AUGUST 1965 


NASA 


GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 


DAMPING OF INTERSTELLAR PLASMA WAVES BY 
THE COSMIC RAY GAS 
by 

Derek A. Tidman* 

Goddard Space Flight Center 
Greenbelt, Maryland 


(Communicated to the Astrophysical Journal) 


(August, 1965) 


^Permanent address: Institute for Fluid Dynamics and Applied Mathematics, University of Maryland, College Park, 
Maryland. 


ABSTRACT 

//&* 

We discuss the propagation of waves in a thermal plasma that co- 
exists with a tenuous gas of relativistic suprathermal particles. Low 
frequency modes such as the magnetosonic wave and Alfven wave ex- 
perience a damping due to the suprathermal particles. Waves propa- 
gating in such regions as the interstellar medium or the plasma of a 
supernova remnant would be significantly damped by the cosmic radia- 
tion fluxes as a result of this mechanism. Although the decay of dis- 
turbances in the interstellar medium (for example produced by stellar 
winds or old supernova remnants) would be hastened by this process, it 
is unlikely that it represents a significant energy transfer to the cosmic 
radiation. For the magnetosonic wave propagating directly across a 
magnetic field the damping has its origin in a cyclotron resonance be- 
tween a wave of frequency u(o > 2 «n o 2 .) and those relativistic particles 
for which fi. = (n oi /y) = where y - (l -v 2 /c 2 )~ 
nonrelativistic ion cyclotron frequency. 
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DAMPING OF INTERSTELLAR PLASMA WAVES BY 
THE COSMIC RAY GAS 

by Derek A. Tidman 

-- : 1 - ■ • ■ * ■■ ' ■ . -r •- , ».•*.? , y 

1. INTRODUCTION t . • 

-■ "V 

Many plasmas in nature contain both a thermal and a suprathermal population of par- 
ticles^ Some obvious examples are the cosmic radiation which coexists with the interstellar 
or interplanetary plasma, or the fluxes of relativistic particles in the Crab nebula. or a 
solar flare, etc. It has been emphasized recently by Parker (1965) that the presence of such 
energetic particles in a thermal plasma can play an important role in the propagation prop- 
erties of waves in the medium. This may be the case even though the number density of 
the energetic particles is much less than that of the thermal plasma since the energy den- 
sity (pressure) of the suprathermal gas may still exceed that of the thermal plasma. k> 

In this paper we point out that a relativistic suprathermal gas also gives rise to an 
important damping of waves in a thermal plasma. This damping is automatically excluded 
if one, uses moment or fluid equations as in Parker's treatment (1965). It emerges natur- 
ally from the relativistic Vlasov- Maxwell equations in a similar way to the conventional 
Landau ^damping (see for example the discussion and references in Chapter 10 of Mont- 
gomery and Tidman, 1964). We have calculated this damping for some "hydro magnetic 
modes," i.e., for waves of frequency much less than the ion cyclotron frequency since such 
waves are thought to be important for energy transport through the interstellar medium. 
However it should be born in mind that such damping rn.echanisms.as we are concerned 
with have their origin in wave-particle resonance effects and are likely to occur in one 
form or another for astrophysical plasma waves oyer much wider frequency ranges, They 
give rise to a damping which is distinct from the viscous damping treated earlier (Parker, 
1955; Piddington, 1957). 

Consider first the mechanism for cyclotron damping (Stix, 1962) of a wave of fre- 
quency +oj t and wave number k propagating at an angle 0 across a magnetic field B 0 . As- 
sume the wave has a nonzero component of its electric vector perpendicular to B 0 . Then 
some of the distribution of particles with velocity v will drift through the wave with just 
the right velocity along B 0 to experience the wave electric field at a harmonic, n, of their 
cyclotron frequency. This occurs if (see Figure 1) 



and gives rise to a resonant transfer of energy between the wave and these particles and 
to a net wave damping. We have only considered the ions for simplicity and fi oi = e|B 0 | /m. c , 
and y - (l - v 2 /c 2 )~ 1/2 . Suppose now one takes the nonrelativistic limit c-»®, y- 1 , and 
considers the low-frequency hydromagnetic part of the spectrum, co T 2 « Cl 2 { . It is then 
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clear since necessarily j v (| j < c that only waves propagating in the range 0 < 0 < 0 M 
= cos' 1 (n 0 ./ kc) can experience such cyclotron damping. In this case there is zero damp- 
ing for waves propagating directly across B 0 (0 - tt/2) as is well known. 

Suppose next we return to consider a relativistic situation in which the suprathermal 
tail extends well into the relativistic region. Then some particles in the range |v |( | < c 
can still satisfy the resonance condition if | (o) r -nfi 0i/ /y)/k cos 0) < c. This remains true 
as 6-n/ 2 if we choose w t - nfi 0i /y. It gives rise to an essentially relativistic damping 
of low frequency waves (^ r 2 « fi 0 2 ) by highly relativistic particles which experience cy- 
clotron resonance with the wave at their reduced cyclotron frequency fij = n oi/ /y = w r < fi oi . 
The physical mechanism for this damping is discussed further following its derivation in 
Equation (32). 

In the following sections we shall calculate the damping decrement for the two cases 
of perpendicular (0 = W2) and parallel (0 = 0) propagation, although of course the same 
effect will occur for arbitrary angles 0. In particular we consider the low frequency 
magnetoaonic wave (& - n/2) and Alfven wave (0 = 0), both of which have phase veloci- 
ties cu,/k = v A where V A is the Alfven speed. If the suprathermal particle energy density 
« becomes important the magnetosonic wave dispersion relation becomes modified. For 
the special case of perpendicular propagation it then corresponds to Parker's "supra- 
thermal mode" (Parker, 1965) with a dispersion relation a> t /k z v A (1 + w/n 0 m. v A J ) 1/1 
where m, and n 0 are the ion mass and thermal number density respectively. 

For the case of hydromagnetic waves propagating through the interstellar medium 
the damping which is produced by the cosmic ray gas depends sensitively on the direction 
of propagation 0. The damping lengths are also frequency dependent and range from a 
value much less than a light year up to indefinitely high values as u r - 0 . The phenomenon 
thus provides a mechanism by which energy is continuously transferred from waves to 



Figure 1— Cone of wave vectors for which cyclotron wave 
damping can occur for a nonrelativistic plasma. 
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suprathermal particles in such regions. In this very wide sense it is thus related to the 
Fermi mechanism - just as the attendant plasma wave damping can be regarded as a 
relativistic form of cyclotron Landau damping. However in the classical Fermi mechan- 
ism particles are visualized as being reflected in collisions with nonlinear wavefronts of 
magnetic field and travel freely between such collisions. In our case the suprathermal 
particles are a part of the plasma and are trapped along a magnetic field B 0 together 
with the thermal plasma. A wave then propagates past these particles through the medium 
continuously depositing energy preferentially into the suprathermal particles through the 
cyclotron mechanism. The waves involved are natural modes of the plasma and one does 
not have to construct models for them. 

It should be pointed out of course that the Fermi mechanism is a more appropriate 
model for the stochastic acceleration of particles in a region of violent nonlinear plasma 
turbulence such as a supernova. Our calculations are appropriate to the dissipation of 
small- amplitude waves in the cosmic ray gas such as might be excited by the old decay- 
ing remnants of a supernova or by stellar winds, etc. Such wave damping however rep- 
resents only a small energy input into the cosmic ray gas compared to that available for 
particle acceleration in the early phase of violent turbulence in a supernova. 

2. BASIC EQUATIONS 

The relativistic Vlasov- Maxwell equations are, 


at 


+ V 




ap 


o 


( 1 ) 


± * - 13B 

ax x t c at 


( 2 ) 


8 n _ HE 4 ?t V 1 f, 

aT X B = C It + e a f a Yd P ' 


(3) 


where f a (p) is the momentum distribution function for the a th species and is normalized 
so that the probable number density n a (x, t) of the a th species is 


= -.j 


f a d P 


where n 0 is the average number density which is taken to be the same for both species 
(we shall only consider an electron-ion plasma for which a = e or i). The relativistic 
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particle momentum is giver* by, 




(, vy“ 

/, P 2 \‘ /2 . 


P - % vy , 

y - 

O'?) ' 

\ + m a C V 

(4) 


Note that the distribution function f a (p) contains both the "thermal" and any relativistic 
"suprathermal" populations of particles which may be present. 

Now in the usual way we consider the behaviour of a small amplitude disturbance in a 
homogeneous plasma in which there is a spatially uniform magnetic field B 0 and for which 
the unperturbed distributions f 0a (|p|) are isotropic functions of p. Thus we write 
f a 3 f 0a + f a (1) , E = o + E (1 > , B = B 0 + B (1) , and linearize in the perturbations. The lin- 
earized Equations (1) - (3) can then be solved for an initial value problem by the method 
of Laplace- Fourier transforms. Defining 


8(k, w) 



-^-r e i (“ ,-k ' x > E<‘> (x, t) 
(*0 3 


(5) 


with lm(o*) > o, one finds a relation of the form R • 8 = I 
for the transform of the electric field where I repre- 
sents initial value terms which we will not write out - 
they are listed in Montgomery and Tidman (1964). 

3. PROPAGATION ACROSS THE MAGNETIC FIELD 
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We choose a co-ordinate system as shown in Fig- 
ure 2 with B 0 along the z-axis and consider a dis- 
turbance for which all the k vectors are along the 
x-axis. The matrix relation for S then simplifies to, 
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Figure 2— Co-ordinate system for prop- 
agation of waves across a uniform 
magnetic field. 
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where there are two cyclotron frequencies, 
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0 a 


Eg [Bpl 
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( 10 ) 


( 11 ) 


and the arguments of the Bessel functions j n are k Vi /n a = kpx/m^ n 0o . It is also important 
to notice that the integrals may go through a singularity of (a>- nfi a ) -1 , and they are de- 
fined as written for im(&>) > o. They can be continued to im(oj) < o in the usual way. 

Now the process of solving for (t) involves inverting equation (6) for 8 and then 
the transforms (5). The field E (1 > then has components proportional to e ‘“ nt from the 
zeros of |R| = o and also some contributions from branch cuts and possible singularities 
of I . We shall assume that the distributions are such that asymptotically for large t the 
disturbance in the plasma resolves itself into modes with frequencies given by the "dis- 
persion relation" |R| = o and that other contributions decay more rapidly in time. For 
more discussion of this point see Montgomery and Tidman, 1964. 

We have not listed R zi in (7) - (9) since it corresponds to a pure electromagnetic 
transverse wave propagating across B 0 with its electric field E (1 > along B 0 , and is not of 
interest to us here. The modes of present interest are those for which e^ 1 ) and E y < are 
nonzero and which are described by the dispersion relation, 

R xx Ryy R xy R yx ~ ® ( 11 ) 


Among these "mixed" modes (i.e., they are generally neither purely longitudinal or trans- 
verse) is one which becomes the magnetosonic wave in the low frequency limit. 
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Figure 3— Schematic plot of a momentum distribution function 
with a relativistic suprathermal tail. 


By inspection of (7) - (9) we see that in the nonrelativistic limit c - «, 0 a = n oa and 
the factor (oj- H a0 ) comes outside the momentum integrations. In this limit there is no 
Landau damping of these waves propagating across B 0 . The damping has its origin in a 
contribution from the singularity of (o>- nO a ) _1 occurring inside the range of momentum 
integration which can only occur in the relativistic case. 


Magnetosonic Waves 

In this section we calculate the damping decrement for magnetosonic waves propa- 
gating through a plasma containing relativistic suprathermal particles. We represent the 
plasma by the distributions, 

s(pi) 

f 0a = K S ( P ll) + C 1 “*> g a (IPI) ^ 1 2 ) 

where the constant * = l ando < (1 -*) « l. The s-function represents a zero- temperature 
’’thermal” plasma containing most of the particles, and the function g a represents the 
tenuous relativistic suprathermal tail (see Figure 3). We have defined 



The constant * has been chosen to be the same for electrons and protons although gen- 
erally g e / g. . This restriction is not important. Further all distributions are normal- 
ized to unity, i.e., 


1 


dpf„a 



1 • 
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Making use of (12) in (7) - (9) the matrix elements become. 



2 TTCO ( 1 “ K) 



Pi 2 ^n*n a 3 Jn ^ 

k 2 v l 2 (a>-nfi o ) 


(13) 


®yx 




Pi 2 <t>x nO a 2 J n J n ' 

kv x (w-nO a ) 


£?a 

^El,(14) 



Now we seek roots a = w t (k) + ia>. (k) of the dispersion relation r xx R yy - R xy R yx = o. 
In order to make this practical we shall make use of three small parameters: 



« l 


(damping decrement is assumed small) 


(16) 



<< 1 (hydromagnet ic wave limit) 


(17) 


(!-/<)« 1 ( suprathermal gas density is small) 


(18) 


Thus making use of (17) and noting that £ " a 2 / fi a0 = 0 for an electron- ion plasma, it fol- 
lows that R xy has two terms of order td 3 /fi a 3 0 and ( l - k ) respectively. Thus the term 
R xy R yx in the dispersion relation contains only terms o|o) 6 /n a ‘ , (i-«) 2 , ( l - * ) o> 3 /fi a 3 0 ] all 
of which we shall neglect compared to the leading terms of r xx R yy . One can verify that 
the roots are thus given approximately by R xx = R yy - o . It should be noted that in ap- 
plying condition (18) we shall also neglect the energy density of the suprathermal particle 
gas compared with b o 2 /s^ although we are aware that the cosmic ray energy density is of 
the same order as the field energy density for the interstellar medium. We are princi- 
pally concerned with a calculation of the damping decrements and allowing for the supra- 
thermal energy density within the framework of the Vlasov equations generates dispersion 
relations considerably more complicated than those obtained by Parker (1965) without 
much alteration of the damping results for the modes we consider. 
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The magnetosonic mode corresponds to the dispersion relation 


and has a wave electric vector along the y axis, and as one easily verifies from Maxwell's 
equations, a wave magnetic vector along the z axis (see Figure lb). The magnetic lines 
of force thus crowd together without bending at the wave crests perpendicular to the di- 
rection of propagation. 

Now the dispersion relation R yy = 0 has the form (noting (17)) 


- , ,2 


"•I® 


+ c 2 k 2 + <v( 1 - k) F(oj) - 0 


(19) 


where 


F (*0 



a 





Pi d Pi n a (Jn ') 2 a Pl 
(a.-nn a ) 


F + i F. 


(20) 


Thus settings = a> r + i <u. in (19) and (20) and equating the real and imaginary parts of (19) 
to zero and dropping terms of o(co. 2 ) or o[a>. (l - *)] we readily find, 



± v. 


yW 0 (rn e "+ m'j) 


(21) 


and for the damping decrement, 


O - x) Fi (*>' + i0) 


(1-k)V a 2 

2c 2 


F ; 


K + i0 ) • 


( 22 ) 


The approximations in the last terms of (21) and (22) follow if the Alfven speed v A satisfies 
v A 2 « c 2 , which it usually does. 

In order to further simplify the expression for <*>. we need the function F i K + i0 )- 
Consider the quantity in (20). This contributes a pair of branch points at 
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Pn = ± i(m a 2 c 2 +p i 2 ) 1/2 which can be connected by a cut in the p N plane passing through the 
point at infinity, and a pair of simple poles at the zeros of (w - r\Ci a ) Regarding Pi as 
fixed for the time being suppose Ph r (p x , n, a) is a value of P[| for which u t -nn a = o . Then 
making a Taylor expansion for p (| = p R we find 


, n v - , . nn 0aP||* 

K" nfi a) - (P|| -P|| ) ~i a 2 ~2 


(23) 


where y R - y(p x , p R ) = nn oa/ ^ r . Thus the Plemelj formulas become in this case, 


' 2 m 2—2 


«."o| 7 1 , - - E . y .V y * m a C (OA) 

1 - FT -nfi a ) - W1 Z- s ( p ir p ii ) Inn p R | l i ne a)HPi) (24) 

— Oa “|| 


where for the sake of definiteness we have assumed^ > o, and l(ne a ) = o if ne a < o or 
= l if ne a > o. The summation £ goes over the two roots p R of « r - nO a = o, which are of 
equal magnitude and opposite sign, and P is the principal value operation. It should also 
be noted that a zero of &> r - nCi a only exists for a finite range o < p x < p 0 which is accounted 
for in (24) by defining j = l in this range and j = o for p x >p 0 where y K = (l +p o ym o 2 c 2 ) 1/2 . 

Now the complex part of F r + iF i derives from the second term of (24). If we note that 
at p 1( = Pl R , r R = nn 0a /u r it follows from (22), (20), and (24) that 


( 1 -K) 77 2 V 


- ^ "l> a ! l fl oJ 2^ I" 1 ! ( ne “) E\ ^ ' f 

a n R ■'-oo Jc 


p opi 2 d Pi ag a 

|p ( J n') 2 5p7 s (prP|i R ) 


(25) 


Next, since p 0 = ( Px 2 + ( P|| R ) 2 y /2 we have 




r*n„ 



(26) 


so that 


Po 



(27) 
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and the quantity d%Jd? L in (25) can be written as (pi/p 0 )(dg o (p 0 )/dp 0 ) because of the 
S -function in the integrand. Thus (25) reduces to 


CO. 
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2(1 -k)7T j V a 2 



+ 00 



I ( ne .) 


g«* (Po) 
Po 



( 28 ) 


FinaUy we note that m e 2 ^ |n 0e | - n^ 2 ^. 2 |nj, j. n (-z) - j n (z), and p 0 = no |n 0a m a |/« r is 

independent of the species a. Thus, 


CO. 
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- 2tt 2 ( 1 


->&)Ttey 


l^oJ 'p 


( 29 ) 


where the dimensionless number ^ depends on the details of the relativistic suprathermal 
particle distributions and is given by, 




® ^ /.i 

2_, I n | [poV (Po) + Po 4 8i' (Po)] 



( 30 ) 


with P 0 - nc |n 0e 


Approximate Expressions and Physical Mechanism for the Damping of 
Magnetosonic Waves 

Consider the last integral in (30). If we are only interested in an approximate result 
we can split the range of integration up as follows, 



( 31 ) 


and use the asymptotic forms of the Bessel function in each range, i.e., 


Jn( 2 > = 



M » l. 


| z | >> n . 


and 


Jn (O 



1*1 « 1 . 


n > 0 - 
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The contribution from the first range is then negligible, and that from the second part of 
(31) gives 



If further, g e ' and g.' are monotonically decreasing functions of p 0 , then due to the 
factor rf 2 in (30) we can approximately neglect all terms except n = l. Thus <u. reduces 
to 


1 d-'ofir) (^-) ln 0 .l p 0 4 [ee'(Po) + gi'(P 0 )] (32) 


where 


Po 



(33) 


The physical mechanism giving rise to this damping decrement can be understood as 
follows. Consider a wave E sino> r t of frequency propagating through the plasma and 
consider only its interaction with the ions for the present. An ion of momentum P 0 is in 
exact cyclotron resonance with the frequency^ = n o . / y 0 , and those particles which are 
nearly in resonance (see Figure 4) have frequencies, 

n oi 

— s a> r ± deo f 

Those particles of slightly less momentum P 0 - AP have a slightly higher cyclotron fre- 
quency £U r If their phase in the electric field E sina> r t is such that they lose energy 

initially they further increase their frequency n oi /y and get rapidly further from exact 
resonance. On the other hand particles with momentum P 0 - £P and initial phase such that 
they gain energy undergo a decrease in n oi /y = which brings them into closer 

resonance for further acceleration. Thus the net result for particles in the range 
P 0 - £P < p < P 0 is a gain of energy from the wave. Similarly the net result for particles in 
the range P 0 <p <P 0 +AP is a loss of energy to the wave. There are more particles in the 
first range than the second which gives therefore a net damping of the wave with a damp- 
ing decrement a> i proportional to the derivative of the distribution function f oi . 
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Figure 4— Schematic plots of the variation of the ion momentum 
distribution and cyclotron frequency. 


4. PROPAGATION PARALLEL TO THE MAGNETIC FIELD - ALFVEN WAVES 

The dispersion relation for circularly polarized cyclotron waves propagating through 
a relativistic Vlasov plasma is well known (see for example Eq. (10.75), Montgomery and 
Tidman, 1964), 


Cl) 


2 



a , 5f 0a (! P I) 
Pi dP + dpj 

y[ w "kv„ ±n a ) 


o 


(34) 


The ± signs refer to right and left polarizations respectively. In the limit co 2 «ft 0 * this 
reduces to the case of Alfven waves propagating along the field B 0 . 

Substituting the distribution functions (12) equation (34) becomes, 



+ 



2 j If* 
pi 2 dpi 

•y (a) - kv h ±n o ) 


= 0 


(35) 
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We shall now make the same approximations as those in (16) - (18) in order to calculate 
the small damping decrement for low-frequency Alfven waves by a tenuous suprathermal 
particle flux. We also assume that the zeros of (oj - kv |( ± fi a ) in (35) occur at values of v N 
such that ( v ,| R ) 2 » - V A 2 , and therefore neglect the real part co t of w in this factor. 

Writing a> = u r + i<a. and dropping small terms (35) then reads, 



Making use of the Plemelj formulas and the fact that g a is isotropic and v A 2 « c 2 leads di- 
rectly to, 


CO 

T 

“k 


+ v 4 


(36) 


and 


- 7T 2 (1 -K) 





(37) 


where the dimensionless number O is given by 


$ 



p 0 p d p [g e (p) 


+ Ki (P)] • 


(38) 


with 


m. fl n . 

i Oi 



(39) 


5. DISCUSSION 

Equations (32) and (37) represent the damping decrements for hydromagnetic waves 
propagating in directions perpendicular and parallel to the magnetic field B 0 respectively. 
It is clear by inspection of these expressions that co. 11 » w. 1 . This is because for parallel 
propagation waves interact with a lower energy part of the suprathermal particle spec- 
trum than for perpendicular propagation. The lower energy particles both gain energy 
more rapidly and are likely to be more numerous. 
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In the following discussion we estimate some characteristic damping lengths for 
waves propagating through the interstellar medium. For this purpose let us consider a 
general integral power law momentum spectrum for the flux of suprathermal protons, 


N(>p) 



sterad 1 


(40) 


We shall assume this is valid for p > p min , and neglect electron fluxes for the present. 
Then equating dN/dp to n 0 v(l - «) g 4 (p) p 2 where v = p/m i y is the speed of a proton of 
momentum p, we have for the momentum distribution function as defined in (12), 


(i - K ) gi (p) 



(41) 


If we define characteristic damping lengths Li and L ufor hydromagnetic waves propa- 
gating perpendicular and parallel to the magnetic field respectively, it readily follows 
from (32), (37), and (41), that 


2 c / c_\ 2 (%A 2 ro M 1 

w |n |V V a) \"e/ ' C ' r ( 3+r )\ W r ) 

1 0 e ' 

provided 



Po 



and 




(42) 


(43) 


(44) 


provided 


P m 



(45) 


The momentum conditions (43) and (45) are simply that those particles responsible for 
the damping should lie in the range p > p min for which (40) is valid. The velocity v 0 is 
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that for an ion of momentum P 0 , i.e., v Q = P 0 /m. > 0 , and v m for an ion of momentum p m 
assuming however that p, - m i v. . In equation (44) if p m is in the relativistic range i.e., 
p a. > m i c > we can replace v m by c to obtain an approximate result for L n . 

Corresponding propagation lengths for waves propagating at arbitrary angles to the 
magnetic field will probably lie between these two values. We see from the conditions 
(43) and (45) that the perpendicular propagation waves are damped only by highly rela- 
tivistic particles, whereas the parallel wave interacts with the low energy part of the 
suprathermal spectrum. We shall next consider as an example the interstellar medium, 
although it should be noted that in regions of stronger magnetic field and more intense 
suprathermal fluxes such as the Crab nebula, or a solar flare, these effects could be 
more important. There are also many forms of such damping for the wider frequency 
range of plasma waves. 

Damping of Waves in the Interstellar Medium by the Cosmic Ray Gas 

For highly relativistic protons the total energy E = pc and equation (40) can be writ- 
ten approximately N - C^m i c 2 /'E) r . The values of r and c appropriate to various ranges 
are given by (Waddington, 1960), 

r 1.4 , C = .5 for 4<pc< 10 4 BeV 
r - 1.5 , C = 1 for 10 4 < pc < 10® BeV 

r = 2 , C = 1,400 for 10® < pc < 10 9 BeV (46) 

Assuming the following parameters for the interstellar medium: n 0 = number density 

for the ionized component = .1 cm -3 , and B 0 = io" 5 gauss, i.e., |fi 0e | = 175, &> e - 1.8 10 4 , 
v A /c = 2.3 10' 4 then leads to 


Li 





light years 


for n 0i >co r > io ' 4 n 0 . . Note that if the suprathermal spectrum parameters F and c re- 
main constant then the numerical factor in this expression for L is proportional to 
(n 0 /B 0 ) . It thus decreases in regions of smaller n 0 or larger B 0 . 

Next consider the case of parallel propagation. Here we have the difficulty that we 
do not know what the very low energy part of the cosmic ray spectrum is in the inter- 
stellar medium. Thus assuming we know the values of C and r for equation (40) only down 
to p min , it follows from (45) that we can only calculate the damping decrement for 


15 


frequencies 



(47) 


Frequencies higher than this would be damped by unknown fluxes of lower energy parti- 
cles with p < p m . . 

* ‘min 

Suppose for example we assume that the values in the top line of (46), i.e., C ^ . 5, 
r = 1.4, are valid down to some value p min . Then it follows from (44) and (45) that Alfvfen 
waves in the frequency range given by (47) are damped in the cosmic ray gas in a length 
scale, 


icr 6 



light years . 


Now we know that the spectrum (40) is a reasonable valid one down to p min c = 1 BeV . 

From this it follows that frequencies in the range w r < 2 . 10' 4 n 0i will be damped in a 
length scale which has a minimum value for oj r = (v A /c)(m ; c/p min ) n oi = 2.10' 4 n 0i of 
L|| = 10" 1 light years. 

Further, suppose we were to assume the spectrum (40) remained valid down to 
Pmin c ~ 100 Mev. This would lead to a damped frequency range o) r < 2 . 10~ 3 fi 0i with a min- 
imum damping length at o> r = 2 . 10' 3 ft 0 . of L y = 10" 3 light years. 

It is thus clear that if there exist any suprathermal particle fluxes in the range 
m i c 2 > pc > a few times the thermal energy, then these would heavily damp a wide fre- 
quency range (cu r <n 0i ) of Alfvfen waves. On the other hand the damping length L x for 
perpendicular propagation involves the highly relativistic particles and is considerably 
longer than L„ , although for the range 10" 2 n 0l < < fi 0 . it is still less than a galactic 

radius. It should be noted of course that the interstellar magnetic field B 0 is nonuniform 
with some characteristic length scale l. Thus strictly speaking our results only apply if 
L,, or L l are much less than i. However even if the reverse situation holds, i.e., I < L (| 
or L L for some frequencies, we would still expect that the same basic process of cyclotron 
damping should occur for waves of wavelength \ << l propagating through such regions. 
These waves would suffer refraction and propagate at varying angles across B 0 . 

Since interstellar plasma waves are damped by the cosmic ray gas it is natural to 
consider whether this contributes a significant energy imput to the cosmic radiation. We 
first note however that the waves discussed here are of small enough amplitude that the 
cosmic ray particles are able to freely traverse the waves with only a minor perturba- 
tion in their trajectories. This is distinct from the reflection of energetic particles by a 
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nonlinear wavefront which is the scattering process basic to the statistical Fermi ac- 
celeration of cosmic ray particles. Suppose further we accept the theory that the bulk of 
the cosmic rays in our galaxy have their origin in the supernova explosions* (Ginzburg 
and Syrovatskii, 1964). It is then clear that there is much more energy available for 
energizing cosmic rays in the early violently nonlinear turbulence described by the 
Fermi mechanism than is available in the final decaying phase when the expanding super- 
nova remnant has degenerated to a region of small-amplitude turbulence. We conclude 
that this wave damping is not likely to be an important source of energy for the cosmic 
ray gas. However it could play an important role in the dissipation of disturbances in 
the interstellar medium caused by old supernova remnants, or stellar winds, etc. 
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* Ginzburg and Syrovatskii (see pages 326-329 of "The Origin of Cosmic Rays") point out that with plausible assumptions for the 
acceleration and escape of suprathermal particles from bounded turbulent plasmas such as supernova shells, one can obtain a 
"universal” integral power law spectrum E -1 * 5 for the energetic particles escaping from such plasmas into the interstellar medium. 
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